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Рассматриваются модели, применяемые для численного исследования процессов с участием минеральной составляющей твердого топлива 
в технологиях энергетического сжигания. 


Введение 

Численное моделирование условий и процессов с 
участием минеральной составляющей при сжигании 
углей в котлах развивается, с одной стороны, как инс- 
трумент оперативного расширения научной основы 
для исследования проблемы шлакования и золового 
загрязнения поверхностей нагрева. С другой стороны, 
возможности моделирования для прогноза последствий 
и корректировки конструкторских и эксплуатационных 
решений стимулируют разработку экспертных систем и 
автоматизированных систем контроля и диагностирования 
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состояния поверхностей нагрева котлов и образования на- 
трубных отложений. Поскольку такие работы невозможны 
без использования математического аппарата, который 
описывает эти сложные и многофакторные процессы, то 
прогресс зависит главным образом от наработок в теории 
физико-математического описания закономерностей тур- 
булентного переноса в двухфазной среде, химического 
реагирования и теплообмена, динамики аэродисперсного 
потока в пучках труб, гетерофазных взаимодействий в 
проницаемом слое и др. Наличие существенного задела 
в области лабораторного и натурного эксперимента по 
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основным процессам горения и аэромеханики способс- 
твовало созданию адекватных математических моделей. 
Анализ состояния дел в этой области [1] констатирует, 
что «математическое моделирование котельного обо- 
рудования становится одним из важнейших способов 
получения наиболее представительной информации об 
аэродинамике, локальном и суммарном теплообмене». На 
этом фоне моделирование процессов с участием минераль- 
ной части (МЧ) выглядит менее проработанным и более 
фрагментарным. Объяснение этому находится, главным 
образом, в том, что эти процессы являются следствием 
всех других многофакторных явлений массо- и теплооб- 
мена при горении и динамике аэродисперсных смесей и 
поэтому наиболее сложны для разработки обобщающих 
математических моделей и доступных для их реализации 
численных алгоритмов. 


Преобразование угольной частицы 

Поскольку одним из основных факторов преобразова- 
ния МЧ топлива является температура, прогнозирование 
свойств минеральных составляющих пылеугольных час- 
тиц возможно с учетом динамики температурных условий. 
Такой подход к исследованию преобразования минераль- 
ных веществ в увязке с температурными изменениями в 
топке и частице разрабатывался В. В. Померанцевым с 
сотрудниками [2,3], причем в одной из работ [3] модели- 
ровались траектории движения частиц, а по зонам столк- 
новения их с ограждающими поверхностями предлагалось 
предсказывать локализацию шлакования или абразивного 
износа. Однако в качестве частиц рассматривались лишь 
обособленные внешние минералы, а их преобразование 
ограничивалось лишь термическим разложением, что 
является начальной стадией в условиях топочного про- 
цесса. На этом, в частности, сказалось отсутствие данных 
о количественных оценках минералогических форм в 
составе МЧ. 

Общее построение модели температурных преобразо- 
ваний МЧ в составе угольной частицы проведено [4] на 
основе теі 6Ѳ^ + 6Ѳ | + б© ' _ 

6Ѳ , 6Ѳ ,6© , 60 , 6© 
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главные составляющие которого известны из работ по 
горению твердого топлива [5]: 
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где: а - коэффициент излучения абсолютно черного 
тела; е - степень черноты; % - коэффициент облу- 
ченности (угловой коэффициент); Т - температур 
облучающей среды. К; Т ч , с ч , Ѵ ч , Ф, р ч - соответственно 
температура, теплоемкость, объем, поверхность, плот- 
ность частицы; Т , Д^ с - соответственно температура, 
коэффициент теплопроводности и теплоемкость газов; 
а\, а',, а' 21 - константы скорости реакций соответс- 
твенно образования СО, СО, и восстановления СО,; 
О,, О,, Оз - тепловые эффекты этих реакций; С я , С 2К 

- концентрации соответственно кислорода и углекис- 

лого газа у поверхности частицы, выраженные через 
концентрацию кислорода; (3! - коэффициент, учитыва- 
ющий изменение числа молей при протекании реакции 
2С + О, = ГГ_ о' , у анта скорости потребления 
кислорода'- ’ 2 7’ 44/32 - мольное соотношение 

СО, и О,. 

Составляющие теплового баланса, характеризующие 
термическое разложение, имеют вид: 

- теплота ж®© с л 5 Оу ч 

т — ( і и I ) I р® ’ ? 

- теплота на® с лб с % л 

і — у ІІІТ I I ) I рі. I , 

где: I - теплота испарения воды; с вп , с л - соответствен)®) 
мкость водяного пара и летучих ингредиентов; 
их плотность; О р;т - теплота разложен)®) 
горючей массы с выделением летучих; с!\Ѵ, сіѴ л - соот- 
ветственно количество влаги и летучих веществ, выде- 
лившихся из частицы за время сіт. 

Теплота, затраченная на преобразования минеральных 
вещест „ 
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а « : п - количество учитываемых компонентов; 
у "оля і-ого компонента в неорганической ма^- 
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теплоемкость минера "ого включения;' 
- тепловой эффект реакции; ,и -доля прореаги- 
ровавшей части і-ого минерального компонента. 

Подстановкой (2-9) в (1) последнее уравнение при- 
водится к виду: 
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Уравнение химической кинетики для реагирующего 
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минерального компонента запишется: 

ІГ 1 — ехр(- р^-)(°\ -ц.)\ 


где п. - порядок реакции, к оі - предэкспонентный мно- 
житель, Е - энергия активации; ц. - текущее содержайй) 
і-ого компонента. 

Уравнение на температуру газовой смеси Т для 
совместного решения с ( 10 ) выводится из теплового 
баланса объема газовой смеси: 
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и представлено [ 6 ] в следующем виде: 
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где С - концентраціи СО в расчетном объеме; 


а 12’ 12 


- соответственно константа скорости и тепловой эффект 
реакции горения СО; 56/32 - мольное соотношение СО 
и СО, к О,; а ф - степень черноты газовой среды; V м 

- относительный выход летучих веществ, отнесенный 
к горючей массе топлива. 

Скорость испарения влаги определено согласно 

_б№а г ( т г - т ( і } 


с* I /А . ВО ч 
8ірі ( + ) 


(12) 


где Т м - температура мокрого термометра; \У - от- 
носительная влажность; А и В - экспериментально 
определенные коэффициенты; 6 0 - исходный размер 
частицы. 

Уравнение для концентрации диоксида угд ода и 
кислорода в расчетном объеме газовой смеси г [ 6 ]: 
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где: к, к 2 - коэффициенты пересчета, учитывающие со- 
ответственно состав летучих веществ и долю диоксида 
углерода от состава летучих веществ в общем объеме 
их продуктов сгорания; С 2 , С - концентрация соответс- 
твенно СО, и О, в расчетном объеме. 

Скорость выделения и сгорания летучих вещфЩ 

горючей массі б с авнениями: 
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Численная модель, представляющая собой с" — ■ 

: і® иТ Т Щ Т ~ А 
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п ’ ’ іеализована [ 6 ] 


вариантов с различными 


,6 


°’ избытками воздуха, 

минеральными включениями. 

Полученные результаты позволяют проводить 
предварительную оценку габаритов камерной топки, 
однако возможности моделирования ограничены из-за 
состояния информации об исходных компонентах МЧ 
в частицах угля. Кроме того, численную модель можно 
сопрягать с аэродинамикой топки только для частиц 
мелких классов, у которых отсутствѵет витание отно- 
сительно потока, что соответствует N 0 = 2 . 


Образование натрубных отложений 

Модель динамики образования золошлаковых от- 
ложений предложена [ 8 ] для условий, где превалируют 
процессы шлакования. За интегральную оценочную 
величину принято используемое в нормативных мето- 
диках значение коэффициента загрязнения: 

І =6 /х , 

представляющего собой термическое сопротивление 
слоя, где: 5 з , - соответственно толщина и коэффици- 

ент теплопроводности натрубных отложений. 

Толщина слоя отложений изменяется от времени и от 
интенсивности про . ия частиц в слое: 

б =| У Л 6 Т, 

(т 2 — т, ) 

где: 2 1 - временной диапазон протекания про- 

цесса, § шл = ОД - интенсивность шлакования, опреде- 
ляемая по формуле [9]: 

г хе, 
л = со і і л, 

где: со ч - [ скорость потока; С ч - концентрация частиц в 
потоке; - коэффициент вероятности столкновений 
частиц с поверхностью; К шл - коэффициент шлакования, 
представляющий собой долю закрепившихся частиц 
относительно столкнувшихся с поверхностью. 

Входящие в (18) величины определятся: 

- скорость золовых и шлаковых частиц в пограничном 
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слое по уравнению [10] 

3® Л® / V 
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00 1 = - 


2*“ (2Х г +Лі )8і 


(18) 


где К - универсальная газовая постоянная, р - дав- 
8 1 ’ молекулярный вес газовой среды, 

’ ' толщина пограничного слоя по пери- 

метру трубы под углом к направлению потока 90° [11]; 
г - радиус трубы; т| [ - кинематическая вязкость газа; Т 

- температура поверхности отложений; 

- концентрация золы в потоке 
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где - расчетный р<; од угля; - зольность р«_ 
чего состояния угля; - объем дымовых газов; 

- доля золы, уносимой из топки в газоходы; 

- коэффициент вероятности столкновения частиц с 
поверхностью трубы [9] 
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где 6 -коэффициент . с , х ;тиц,п- показатель 
1 Д — О / Ои «Е, 

полидисперсности и - относительный 

размер частиц, определяемые по уравнению Розина- 
Раммлера [9]; 

- коэффициент шлакования [8] 
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где - радиационный и конвективный тепловой по- 

ток, воспринятый поверхностью нагрева, р з - плотность 
золового отложения, К - коэффициент, учитывающий 
энергетический баланс взаимодействия с поверхностью 
нагрева золошлаковых частиц как совокупности твердых, 
размягченных и расплавленных частиц. 

Коэффициент К в (20) зависит от разницы темпе- 
ратур между газовым потоком и натрубными отложе- 
ниями, от скорости частиц, от теплопроводности газов 
и натрубных отложений, от реологических свойств. В 
свою очередь, реологические свойства рассматриваются 
как функция минерального состава частицы, состава 
газовой среды и температуры частицы, для расчета 
которых используются уравнения системы ( 1 0 02^) 
Поэтому применение модели осложнено теми же огра- 
ничениями, что были отмечены для данной системы. 
Кроме того, расчет концентрации частиц в потоке по 
(19) дает среднее по топке значение, в основе которого 
лежит нормативная величина выноса золы из топки 
й _ | мая из о" ""ращенного золового баланса 
Здесь представляет собой коэффи- 
циент шлакоулавливания, который зависит от фракци- 
онного состава дисперсной среды, аэродинамической и 
температурной обстановки в зоне активного горения, а 
поэтому связан с режимными факторами. 

Для повышения эффективности моделирования за- 
грязнения поверхностей нагрева следует моделировать 


локальную концентрацию золовых частиц в потоке, а это 
связано с расчетом процесса шлакоулавливания. 

Другим важным ограничением в рассматриваемой 
модели является отсутствие увязки коэффициента 
теплопроводности натрубных отложений с трансфор- 
мацией слоя, выражающейся в изменении состава и 
плотности. 

Параметры пылегазовой среды 

Один из подходов к оценке полей температур, 
скорости, концентрации двухфазного потока в топке 
и газоходах котла представляет собой использование 
численной методики ТГУ [12]. Математическая мо- 
дель состоит из двух систем уравнений. Одна из них 
предназначена для исследования пространственного 
стационарного турбулентного движения газовой среды 
в объеме топки, базируется на уравнениях Рейнольдса 
и - е - модели турбулентности. Способ задания 
граничных условий для системы уравнений опреде- 
ляется типом границы. Для входных границ (сечение 
ввода первичного и вторичного воздуха) используются 
однородные распределения для всех характеристик те- 
чения. На выходе из расчетного объема ставятся мягкие 
граничные условия, причем расположение выходной 
границы должно быть достаточно удалено от зон воз- 
можной рециркуляции потока. В качестве граничных 
условий на стенках топочной камеры используется 
метод пристеночных функций Лаундера-Сполдинга. 

Модель была протестирована на экспериментальных 
данных визуализации течения, полученных в ТПУ при 
гидравлическом моделировании топки котла БКЗ-220- 
1 00 ЖШ . Результаты тестовых расчетов дают основания 
применять данную модель для численного исследования 
аэродинамики топок. 

Другая система уравнений описывает движение 
частиц в рамках лангражева подхода с учетом дина- 
мического теплового скольжения относительно газа 
и с учетом их горения. При этом частицы считаются 
сферическими, а столкновения их между собой не учи- 
тываются. В качестве начальных условий используются 
значения величин, соответствующие местоположению, 
скорости, температуре и размеру частицы в выходном 
сечении горелки. 

Основные системы уравнений могут сопрягаться 
с моделью теплообмена и течения жидкой шлаковой 
пленки, образующейся на стенках камеры сгорания 
топки [13]. С учетом этого дополнения оказалось воз- 
можным исследовать эффективность управляющего 
воздействия как на поведение топливо-воздушной 
смеси, так и на распределение МЧ между шлаком и 
уносом от ряда конструктивно-режимных факторов 
топочного процесса: сочетания включенных в работе 
горелок, ориентации горелочных струй, тонкости по- 
мола, конфигурации камеры. 

Заключение 

Состояние проблемы численного моделирования 
процессов сжигания углей характеризуется наличием 
адекватного физико-математического описания приме- 
нительно к отдельным стадиям поведения минеральной 
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части топлива в котельных агрегатах. Повышение 
эффективности моделирования связано с построением 
укрупненных моделей на базе частных, что оправдано 
после повышения эффективности имеющихся фраг- 
ментов моделирования за счет устранения недостаточ- 
ности или неопределенности информационной базы по 
минеральной части углей. Прежде всего это относится 
к оценке содержания компонентов, термофизическим 
характеристикам и параметрам их изменения. Необходи- 
мость их экспериментального исследования как основ- 
ного пути решения вопроса по-прежнему актуальна. 
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